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Abstract. Cosmic voids have been shown to be an effective probe of cosmology, complementary 
to galaxy clusters. In this work, we present a simple theoretical framework for predicting 
of the size of the largest voids expected within a given redshift and volume. Our model is 
based on the exact extreme-value statistics which has previously been successfully applied 
to massive galaxy clusters. We implement our formalism using the void-abundance models 
and compare the extreme-void predictions to simulations and observations. We find that the 
simplest void models can only explain the extreme-void abundance with ad hoc parameter 
adjustments. We argue that extreme-void distributions should be used as an additional test 
on theories of void abundance. 
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1 Introduction 

In the standard concordance (ACDM) cosmology, cosmic structures on the largest scales 
consist primarily of massive galaxy clusters (typically of masses 10^^Mq) and large cosmic 
voids of radius ~ 100 Mpc. Thanks to ambitious ground and space-based surveys and a 
suite of large V-body simulations, the distribution of massive galaxy clusters have now been 
studied in great detail. Several analytic and semi-analytic mass functions are now available 
to describe and predict cluster distribution [1-5]. These prescriptions have generally been 
shown to be in good agreement with simulations and observations. 

Cosmic voids, on the other hand, have not received the same level of attention as 
clusters. The main difficulty is the lack of consensus of what defines a void, for instance, the 
underdensity threshold below which a region is said to be a void. Furthermore, whilst clusters 
can be observed via the mass-luminosity relation, true cosmic voids are empty of dark matter, 
and so cannot be directly observed. The relation between galaxy voids and dark matter voids 
is still not completely understood (see [6, 7] for recent progress). 

Nevertheless, a number of recent simulations have shown convincingly that voids, too, 
are a sensitive probe of cosmology, because they preserve much of the conditions in the 
early Universe. Using indirect probes such as gravitational leasing and the Integrated Sachs- 
Wolfe effect, their abundance, shape and density profile have been shown to be effective 
probes of dark energy [8-10], dark matter [11, 12], gravity theories [13-15] and primordial 
non-Gaussianity [16, 17]. 

The first convincing analytic model of void abundance was given by Sheth and Van de 
Weygaert (SVdW) [18]. Although it has since undergone minor modifications by a number of 
previous authors [19-22], the heart of the theory, based on excursion-set formalism, remains 
largely unaltered. In this work, we ask: how large are the largest voids predicted by the 
SVdW theory? Our focus on the largest voids is based on the fact that they have the most 
interesting phenomenology, and may be responsible for large-scale anomalies in the CMB 
anisotropies, particularly the notorious Cold Spot [23-26], and may produce similar anomaly 
in the 21cm fluctuations [27]. 
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At the core of our work is the exact extreme-value statistics previously used for predicting 
the distribution of the most massive clusters [28, 29]. We will use this framework, together 
with the SVdW theory, to produce a probability distribution of the largest voids expected 
to be observed at a given redshift. Our framework is then tested against simulation and 
observational data. We then discuss how this approach can be improved and extended to 
other void models. 

Even though simulations give us access to the full distribution of void number counts, 
the results are highly dependent on the choice of void finder {e.g. ZOBOV [30], VIDE [31], 
amongst others, see [32] for a recent discussion). Our goal is in this work is to investigate 
whether the extreme-void predictions from current theory is consistent with the extreme voids 
seen in simulations and observations. And if not, how could we address the inconsistency? 

Throughout this paper, we will use the Planck-l-WMAP cosmological parameters as 
given in [33]. 


2 Void abundance 


We begin with an overview of the SVdW theory, which gives the differential abundance of 
voids in the radius interval [i?, i? -|- dR\. The model is based on the so-called ‘excursion set’ 
theory, in which the local linear overdensity, (5(x), performs a random walk in the (<5, cr^) 
plane, where cr^ is the variance of the overdensity, smoothed on comoving scale R: 

^k^P{k)W^{k,R). ( 2 . 1 ) 

D{z) is the amplitude of the growing mode (normalised to D(0) = 1), P{k) is the linear 
matter power spectrum and W{k, R) is the Fourier-space window function. A void of radius 
R is deemed to form when 5 first makes an ‘excursion’ below the barrier 6^. In the original 
paper, the value 6y = —2.7 was recommended (this corresponds to a nonlinear overdensity 
~ —0.8). This value of 5y is derived from the ‘shell-crossing’ condition, analogous to the 
Einstein-de Sitter critical overdensity, dc = 1.686, commonly used to define an object formed 
by spherical collapse. 

The SVdW model gives the differential abundance of voids with radius R as [18] 


/(cr) diner 


-1 


dn 


din7? Vs{Rl) d\nR l 


Rl=R/1.7 


f{a) = 2 'kx^ je s\u{j'KV), 

V = {1 + 5c/\5y\)~^ , x = 'Da/\5y 


( 2 . 2 ) 

(2.3) 


where Vs{r) = jd. The oscillatory behaviour of / at high R can be avoided with the 
approximation given by Jennings, Li and Hu [19]: 


x< 0.0276, 

y TT (T ’ — ’ 

27rx^ sin(j7rT>) otherwise. 


(2.4) 


The SVdW theory has been shown to be consistent with dark-matter simulations where 
only spherical voids are considered [19]. However, void finders based on the watershed 
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algorithm typically pick out voids that are non-spherical and irregular. This led many authors 
to take 6v as a free parameter in the model, and simply adjust 5v to match simulation results. 
Whilst this is a rather ad hoc approach, such an adjustment appears to do the trick, for 
instance, in accurately matching the void abundances in the simulations of [34], with the 
parameter h.,; ~ — 1 for z <1. We will return to these simulation results in §6. 


3 Void number counts 


The probability density function (pdf) for voids with radius in the interval [R, R + dR], in 
the redshift bin centred z, width Az, is given by 


f{R,z) 


/sky dn 

Ntot Jz-Az/2 dz i?dln/?’ 


(3.1) 


where /sky is the fraction of the sky observed and dV/dz is the volume element given by 


dV_ dvr 
d^ “ 


f 


dz' 


H{z') 

H{z)^Ho[Qm{l + zf + nAY^\ 


(3.2) 

(3.3) 


In (3.1), the total number of voids, Wot, with radius above Rmin in the redshift bin [z — 
Az/2, z + A2;/2] is 


-Wot(^) — /sky / 
JR 


^max,V dR rz+Az/2 

- / (aZ -. 

R Jz-Az/2 dzdlni?' 


(3.4) 


where i?max,v is determined by the volume of the redshift bin. 

Fig. 1 shows total number of voids Wot as a function of z (top panel), with Az = 0.05, 
-Rmin = 1 and /sky = 1- Since we will be concerned with the largest voids, our results are 
insensitive to changes in Rmin- In the same plot, we consider what happens when 6^ varies 
from —0.9 to —1.1, and the fractional difference compared to the 6^ = —1 model is plotted in 
the lower panel. 

We see that when the void barrier is decreased (voids become emptier), the number 
count increases (in this case by about 10 — 15%), and vice versa. This can be interpreted as 
the fragmentation within the void population - smaller voids are typically emptier, and they 
can combine to form shallower, larger voids, which are therefore less numerous. This simple 
picture is consistent with detailed studies of void profiles in previous simulations [7, 13]. 

Having obtained the pdf (3.1), we can construct the cumulative probability distribution 
(cdf), R(R), which gives the probability that an observed void has radius < R. It is given by: 


R(R) = 



/(r) dr. 


(3.5) 


4 Distributiou of extreme voids 

Harrison and Coles [35] showed that the cdf naturally leads to the prediction of the most 
massive clusters expected at a given redshift. This section follows their treatment, but adapted 
to cosmic voids. 
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Figure 1. The top panel shows total number of voids fVtot and its variation with respect to the void 
underdensity Sy. The lower panel shows the fractional difference when compared with the fiducial 
model with Sy = —1. 


Consider N observations of voids drawn from the cdf, F(R). We ask: what is the 
probability that the largest void observed is of radius R*7 The cdf required is 

<^{R\N) = Fi{R<R*)... Fn{R < R*) = F^{R*), (4.1) 

where the last step assumed that void radii are independent, identically distributed variables 
(we will come back to this point). Differentiating <1> with respect to R* gives the exact 
extreme-value pdf 

(^{R\N) = = Nf{R*)[F{R*)f-\ (4.2) 

Figure 2 shows an example of the exact extreme-value pdf, 4>{R), at z = 0.225. The pdf 
is non-Gaussian but can be well-approximated as a sum of a few Gaussian terms. On the 
figure we also show the contour in dashed lines representing the 5th and 95th percentile of 
the distribution. 

Assume that all available voids are observed (N = Atot)) the peak of the pdf occurs where 

^<(-(ii*,Wot) = 0. (4.3) 

It follows that the peak of the extreme-value pdf occur at the zero of the function 

X{R) = iN-l)f + F%. 
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Figure 2. The extreme-value pdf, (j), a.t z = 0.225. The values of the void radius, i?, corresponding to 
the 5th (left) and 95th (right) percentiles are shown in dashed lines. 



Figure 3. The distribution of the peak of the pdf, (/^(i^max): where i?max is the largest void expected, 
plotted as a function of redshift. We also show its variation with the void underdensity 6y. Note that 
the ‘shell-crossing’ value (5„ = —2.7 predicts an unfeasibly small radius of extreme voids. 

These zeroes can be used to roughly locate the most likely extreme radius (although when 
comparing with data in the next section, we will again use the percentiles). 

Fig. 3 shows Rmux, the peak of the extreme-value distribution (i.e. the root of X{R) = 0), 
as a function of redshift. With 6^ = —1, we see that the largest void predicted by the SVdW 
theory is most likely ~ 50 /i“^Mpc in radius, at z ~ 0.2 — 0.3, with i?max very gradually 
decreasing at higher redshifts. This slow redshift evolution agrees with the simulation of [36], 
in which the maximum void radius appears to be almost independent of redshift. 

4.1 Sensitivity to Sy 

In Fig 3 , we also plotted the extreme-void predictions for varying Sy. We observe that 
decreasing results in larger extreme-voids across all redshifts, in accordance with the 
previous section where we deduced that larger voids are shallower. 

In addition, we also included the case 6y = —2.7, the shell-crossing value originally 
recommended by SVdW. However, such a large magnitude of predicts the peak of i?max 
at only ~ 20/i“^Mpc, which is effectively ruled out by recent simulations in which much larger 
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Figure 4. The sensitivity of i?max(z) to dark energy. Left: The dark-energy density parameter Oa 
is varied within 10% of the fiducial value (whilst flatness is maintained). Right: The dark-energy 
equation-of-state parameter w (assumed redshift independent) is varied from —0.8 to —1. 


voids have been seen (we will elaborate on this in the next section). Our finding adds weight 
to the poor fit of (5^, = —2.7 to the void abundance dn/dlnR as seen in the simulations of [34]. 

4.2 Sensitivity to dark energy 

The radius of extreme voids is linked to the properties of dark energy. Here, we consider the 
effects of varying the dark-energy density parameter, Oa, and equation of state, w, from the 
fiducial values (Oa = 0.73, u) = —1). 

If Oa were to increase whilst flatness is maintained, then the matter density parameter 
must decrease. This gives rise to a reduced ag and, intuitively, we expect structures to be 
generally less pronounced; clusters are less massive, and voids are smaller. 

If the dark-energy equation of state, w, is larger than —1, we expect that void radius 
will not grow as fast at late times, in accordance with Friedmann equation which predicts 
7?~t2/(3(v;+i))_ Hence, we expect smaller voids if rc > — 1. 

Figure 4 shows the peak of the extreme-void pdf when and w are changed from their 
fiducial values (Ha is varied whilst flatness is maintained). These behaviours are in accord 
with our intuitive discussion above. Note that the effect of varying w seen here qualitatively 
agrees with the results in [37] , in which the effects of w on voids are further investigated. 

In both cases, the variations shown are degenerate with small changes in 6^ in the fiducial 
model. However, current observational constraints on Ha and w are so tight that we do not 
expect these degeneracies to be significant. 

5 Comparison with simulations 

We now compare our extreme-void model with simulation results. We study the void 
distributions in two simulations: 1) Chan, Hamaus and Desjacques [34], using a 1.5 /i“^Gpc 
box, containing 1024^ particles, at redshift snapshots z = 0,0.5 and 1, 2) Sutter et al, using 
a 1 /i“^Gpc box, containing 1024^ particles, at snapshot z = 0. The maximum void were 
extracted from these distributions, giving 4 data point in the z — TZmax plane. 

On top of these data points, we superimposed the theoretical extreme-value distribution, 
calculated as in §4, except that the pdf and the total number of voids are calculated without 
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Figure 5. The radius, i?max, of the largest void computed using the SVdW theory without volume 
factors. The diagonal stripe comprises the peak of the extreme-value pdf (centre - solid line), and the 
95% and 5% percentiles of the distribution (top and bottom dashed lines). Data points (x) are from 
the simulation of [34] whilst the point (*) is from Sutter et aZ.[39]. 


redshift binning: 

(5.1) 

(5.2) 

Fig. 5 shows the result with 5.^ = —1 in solid line, with dashed lines representing the 
5th and 95th percentiles of the distribution. 

We see that the data from Sutter et al. fall within the central 90% of the distribution. 
However, those from Chan et al. fall significantly outside the contours. Informed by the 
results in Fig. 3, it appears that their simulation prefer (!.„ > — 1 at z ~ 0, and decreasing 
with redshift. 

This discrepancy from their recommended value of ~ —1 can be explained by examining 
their Fig. 2 carefully. Even though the bulk of the void distribution can be described by the 
SVdW model with 6v ~ —1, the mismatch at high values of R appears to be significant. 

In conclusion, our results favour a redshift-dependent barrier <5^(2;) - similar to the 
modifications of the SVdW model advocated by [21, 38]. Although SVdW theory with 
(5^ ~ —1 appears to fit dn/dhiR well across the redshifts, we have shown that the data on 
extreme voids are not consistent (be. outside the central 90%) with theoretical expectations 
at any of the redshift snapshots. This illustrates the use of our extreme-value formalism as 
an additional consistency check for void theories. 

The reader may be aware of much larger void radii that have been reported in other 
simulations and observations. For instance, [26, 36] reported voids as large as 350 h~^ Mpc, 
which, according to the SVdW, should be unobservably rare. However, these voids are not 
dehned using the same criterion as the SVdW definition (in which voids are region with mean 
density pv < 0.2p). It is beyond the scope of this work to try to reconcile these two definitions. 
We simply note here that such enormous voids cannot yet be theoretically described using 
the SVdW framework without further ad hoc adjustments. 


f{R,z) = 


I dn 




tot 


poo 

{z)= / 


tot ' 
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Figure 6. The radius, i?max, of the largest void expected in each redshift bin with Ax = 0.05 
computed using the SVdW theory . The diagonal stripe comprises the peak of the extreme-value pdf 
(centre - solid line), and the 95% and 5% percentiles of the distribution (top and bottom dashed lines). 
Data points (x) are from the SDSSIO catalogue [40]. 

6 Comparison with observation 

We perform a similar comparison with voids identified from the SDSSIO catalogue as presented 
in [40]. Voids in the catalogue were binned with Ax = 0.05, and the maximum in each bin 
is picked out. Changing the bin size has no significant effect on neither the theory nor the 
observed extreme-void distributions. The result is shown in Fig 6. The theory curve with 
5v = —1 (solid line) is calculated using the formalism in §4. The 5th and 95th percentiles are 
shown in dashed lines. 

As in the previous section, it appears that no single value of 6y is a suitable fit to the 
observed extreme voids. In contrast with the simulated voids, real voids appear to favour 6{z) 
which increases with redshift. 

The discrepancy between theory and observation is not surprising: it is unclear how the 
SVdW predictions, which are relevant only for dark matter voids, can be extrapolated to 
galaxy voids. Furlanetto and Piran [20] previously tackled this question by modifying the 
SVdW theory to fit voids in the 2dF survey. Perhaps a formalism such as theirs is could be 
used our extreme-value formalism to ameliorate the discrepancy seen in Fig. 6. 

7 Conclusion and discussion 

In summary, we have given a simple theoretical framework for the prediction of the size 
of the largest voids expected within a given redshift bin. Our model is based on the exact 
extreme-value calculations previously applied to galaxy clusters. We implement our formalism 
using the void-abundance model of Sheth and Van de Weygaert, with 5y taken as a free 
parameter. Our results show that the simplest SVdW with 5v = a constant is inconsistent 
with the extreme sizes of voids seen in the simulations and observations. For instance, 6y 
which decreases with redshift would produce a better fit to simulations. 
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Figure 7. Comparison between the SVdW model and the volume-conserving ‘Fdn’ model of [19]. 

The central message in this work is that any theory which predicts void abundance must 
also give a sensible prediction for the distribution of extreme voids. Our work provides a 
framework to perform this check. 

There are several scopes for improvement of our extreme-value model. 

• Void correlation. Of course, voids are not completely uncorrelated, as previously shown 
in the theoretical work of [21, 22, 34, 41]. Preliminary works suggest that introducing 
a void correlation will reduce the SVdW abundance for spherical voids, although the 
effects on the largest voids are likely to be small and predictable by linear theory [42, 43]. 
Our work on including void correlation into the extreme-void statistics is in progress. 

• Alternative void abundances. There is increasing evidence that the SVdW model cannot 
accurately explain simulation data without ad hoc adjustments, simply because real 
voids are unlikely to be spherical, nor are they all likely to form out of the shell¬ 
crossing condition. Apart from simple modification of the SVdW model, alternative void 
abundances include those directly derived from the Press-Schechter cluster abundance 
[16, 44], or one derived using de-Sitter geometry [45]. The extreme-value formalism 
presented here remains unchanged and can be readily applied using these alternative 
models. 

A crucial shortfall of the SVdW model is the fact that it assumes that voids are 
isolated spheres, and as voids evolve from the initial Lagrangian space to the late-time 
Eulerian space, their volume is not conserved (instead, their comoving number density 
is conserved - and only their sizes change). Jennings, Li and Hu [19] showed that 
this assumption can yield volume fraction in voids exceeding unity. They suggested a 
modified volume-conserving SVdW model - so-called ‘V dn’ model, and showed that 
this formalism produces a good fit to their simulation, but only for spherical voids with 
R < 2Qh~^ Mpc (using (5^=-2.7). 

It is curious to see if the V dn model performs any better than the SVdW model when 
incorporated into our extreme-value formalism. Figure 7 shows the comparison (with 
(5„=-l in both cases). Unfortunately the V dn suppresses the radius of extreme-voids, 
and exascerbates the underprediction of extreme void radii. 

• Alternative void definitions. Voids are defined in the SVdW theory as spherical regions 
in which the density satisfies py < 0.2p (where p is the mean matter density). We 
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have found that decreasing the ratio Pvl~P to below 0.2 increases the extreme sizes of 
voids. However, such an adjustment is unsatisfactory as there are consistency relations 
between this ratio and 5.^ (see the Appendix of e. 5 . [19]). Even with the consistency 
relation in place, we were unable to reproduce a significant number of extreme voids 
with radius beyond ~ 100 Mpc unless model parameters describe extremely shallow 
voids. We anticipate that this problem requires a fresh approach to the modelling of 
void abundance, which must include the shape and profile of voids. 

Perhaps the greatest challenge for void theory is the prediction for the abundances of real 
galaxy voids seen, for instance, in the SDSS data [46, 47]. A solution to this problem will have 
a significant implication, in particular, on explaining large-scale observables associated with 
supervoids. Given the highly nonlinear and complex nature of real voids, perhaps it is more 
productive to strive towards a semi-analytic modelling of voids led by accurate simulations, 
similar to what has been achieved with the universal-mass-function approach for galaxy 
clusters. 
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